Complete measurement of three-body photodisintegration of 3 He for photon energies 

between 0.35 and 1.55 GeV 
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The three-body photodisintegration of 3 He has been measured with the CLAS detector at Jefferson 
Lab, using tagged photons of energies between 0.35 GeV and 1.55 GeV. The large acceptance of the 
spectrometer allowed us for the first time to cover a wide momentum and angular range for the two 
outgoing protons. Three kinematic regions dominated by either two- or three-body contributions 
have been distinguished and analyzed. The measured cross sections have been compared with results 
of a theoretical model, which, in certain kinematic ranges, have been found to be in reasonable 
agreement with the data. 

PACS numbers: 21.45.+V, 25.20.-x 



I. INTRODUCTION 

The study of the electromagnetic properties of the 3 He 
nucleus is the optimal starting point to assess the im- 
portance of many-body interactions between nucleons in 
nuclei [1, 2]. In particular, in the 7 3 He — > ppn reac- 
tion, threc-nucleon currents dominate in certain regions 
of phase space [3-5]. In fact, a pp pair has no dipole 
moment with which to couple and the charge- exchange 
current vanishes within a pp pair, so that the one- and 
two-nucleon currents are suppressed in those regions. 
The small number of nucleons involved makes possible 
kinematically complete experiments, and exact Faddeev 
ground-state wave functions, as well as exact wave func- 
tions for the continuum three-body final state at low en- 
ergies (below the pion-production threshold), are avail- 
able [6, 7]. 

Although the calculations of the 3 He ground-state 
wave function have reached a high level of accuracy in 
reproducing the bound-state properties [6, 7], the calcu- 
lation of the continuum three-nucleon wave function is 
less developed at higher energy; a full treatment of the 
three-body photodisintegration of 3 He has been possible 
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only at energies E 7 < 300 McV. As the energy increases, 
the number of partial waves and open channels becomes 
very large and, so far, no calculations that are both ex- 
act and complete have been done in the GeV region. Not 
only would a very large computational effort be required 
to do so, but also a treatment of the absorptive part of 
the nucleon-nucleon interaction (coupling to other open 
channels that is not taken into account in potential-based 
calculations) should be implemented. 

A different approach has been taken by Laget [3, 8-11], 
who has employed a diagrammatic model for the evalu- 
ation of the contribution of one-, two-, and three-body 
mechanisms in the cross section for the photodisintegra- 
tion of 3 He. Rather than relying on a partial-wave expan- 
sion, this approach relies on the evaluation of the dom- 
inant graphs whose amplitudes are related to one- and 
two-body elementary amplitudes. The parametrization 
of these elementary amplitudes incorporates absorptive 
effects due to the coupling with other channels, which be- 
come more and more important as the energy increases. 
The comparison of these model predictions with exper- 
imental data provides us with a good starting point to 
understand the nature of three-body interactions in 3 He 
for photon energies in the GeV region. 

At stake is the link with three-body forces. In the 3 He 
ground state, three-body forces involve the exchange of 
virtual mesons between nucleons and the creation of vir- 
tual baryonic resonances. The incoming photon can cou- 
ple to each of these charged particles. Below the pion- 
photoproduction threshold, all the particles remain vir- 
tual and the corresponding three-body meson-exchange 
currents (MEC) contribute only weakly to the cross sec- 
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tion. When the photon energy increases above the vari- 
ous meson- or resonance-production thresholds, these vir- 
tual particles can become real - they can propagate on- 
shell [1]. The corresponding sequential scattering ampli- 
tudes are considerably enhanced and can dominate cer- 
tain well defined parts of the phase space. Kinematically 
complete experiments allow one to isolate each of the 
dominant sequential rescattcring amplitudes. They ana- 
lytically reduce to three-body MEC at lower energy, and 
put constraints on the corresponding three-body current. 

Several low-energy (<100 MeV) experiments have been 
performed since the publication of the results of the first 
measurement of the three-body photodisintegration of 
3 He in 1964 [12]. But only a few have been performed at 
intermediate photon energies up to 800 MeV, in limited 
kinematics [13-15] as well as with large-acceptance de- 
tectors [16-18]. They show good agreement with Laget's 
predictions provided that the 37V mechanisms, based on 
sequential pion exchanges and A-resonance formation, 
are included in the calculations. Since these mechanisms 
dominate well defined parts of the phase space, a better 
understanding of the nature of many-body interactions 
requires one to perform a high-statistics full Att inves- 
tigation, probing the three-body breakup process for all 
angular and energy correlations of the three outgoing nu- 
cleons. Also, the extension to the high-energy (75 7 > 1 
GeV) region, where no experiment has been performed 
until now, can be expected to open a window on other 
kinds of many-body processes. 

This paper reports on a measurement of the three-body 
photodisintegration of 3 Hc performed in Hall B at Jeffer- 
son Lab [19]. Photon energies between 0.35 GeV and 1.55 
GeV were used, and wide angular and momentum ranges 
for the outgoing particles were covered. These features, 
along with the high statistics collected, allow us to se- 
lect the most interesting two- and three-body processes, 
to compare their relative importance, and to determine 
their variation with photon energy. 

The experimental setup is described briefly in Sec. II, 
the salient points of the data analysis in Sec. Ill, and 
the model calculation in Sec. IV. Our results for several 
kinematic regions are presented in detail and compared 
with the model calculation in Sec. V, and summarized in 
Sec. VI. 



II. EXPERIMENTAL SETUP 

The experiment was performed at the Thomas Jeffer- 
son National Accelerator Facility, in Hall B, using the 
CEBAF Large Acceptance Spectrometer (CLAS) [20] 
and the bremsstrahlung photon tagger [21]. The electron 
beam energy was 1.645 GeV, corresponding to two passes 
of the CEBAF accelerator; the current was 10 nA dur- 
ing regular production runs and 0.1 nA during tagging- 
efficiency calibration runs. The photon beam was pro- 
duced by the electron beam striking the radiator, a thin 
layer (~ 5 x 10 -5 radiation length) of gold deposited on a 



thin carbon backing, which was placed 50 cm before the 
entrance of the tagger magnet. The electrons interact- 
ing in the radiator were deflected by the magnetic field 
of the tagging magnet, and those with energy between 
20% and 95% of the incident electron beam energy were 
detected by two layers of scintillators (E-counters, mea- 
suring the energy of the electron, and T-counters, mea- 
suring its time [21]) placed in the magnet focal plane. 
Thus, photons in the energy range from 0.35 to 1.55 GeV 
were tagged. Two collimators were placed in the beam- 
line between the tagger and the 3 He target, in order to 
eliminate the tails from the photon beam and to give a 
small and well defined beam spot on the target. The 
data were obtained using a cylindrical cryogenic target, 
18 cm long and 4 cm in diameter, filled with liquid 3 He 
and positioned approximately 20 m downstream of the 
tagger radiator in the center of the CLAS. A lead-glass 
total absorption counter (TAC), almost 100% efficient, 
placed approximately 20 m downstream from the center 
of the CLAS detector, measured the tagging efficiency 
during low-flux calibration runs. 

The CLAS is a magnetic toroidal spectrometer in 
which the magnetic field is generated by six supercon- 
ducting coils. The six azimuthal sectors are individu- 
ally equipped with drift chambers for track reconstruc- 
tion, scintillation counters for time-of-flight measure- 
ment, Cerenkov counters for electron-pion discrimina- 
tion, and electromagnetic calorimeters to identify elec- 
trons and neutrals. The polarization of the CLAS torus 
was set to bend the negatively charged particles toward 
the beam line. In order to achieve a good compromise 
between momentum resolution and negative-particle ac- 
ceptance (required by other simultaneous experiments) 
the magnetic field of the CLAS was set to slightly less 
than half of its maximum value, corresponding to a torus 
current of 1920 A. A coincidence between the tagger and 
the time-of-flight scintillators defined the Level- 1 trig- 
ger for accepting the hadronic events. For the first time 
in CLAS, a Level-2 trigger, which selected the events 
from among those passed through Level- 1 that have at 
least one "likely track" in the drift chambers, was also 
used [20]. More than a billion events of production data 
were obtained with 3 He (plus a few million events taken 
with the target empty), at a data- acquisition rate slightly 
greater than 3 kHz. 



III. DATA REDUCTION AND ANALYSIS 

A. Channel Identification 

In order to isolate the ppn channel, a pp coincidence 
(with no other charged particles) in a time window of 
±1 ns with a tagged photon defined the minimum con- 
dition for an accepted event, since the time interval be- 
tween beam pulses is 2 ns. This coincidence time is shown 
in Fig. 1 for a subset of the raw data. The two protons 
were identified by their mass, deduced from their momcn- 
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FIG. 1: Coincidence time for a subset of the raw data. The 
vertical lines in the upper panel indicate the time window 
for accepted events. Random coincidences from neighboring 
beam pulses are visible in the lower panel. 



turn measured in the drift chambers and their velocity 
measured with the timc-of- flight scintillators, as shown 
in Fig. 2. 
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FIG. 2: Velocity j3 — v/c spectrum, as a function of particle 
momentum, for charged particles detected in the CLAS. 

A cut on the interaction vertex, based on the analy- 
sis of empty-target runs, was performed to eliminate the 
background from pp events originating outside the target 
volume. Eliminating the events having the z component 
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FIG. 3: Distribution of the z component (along the beam 
line) of the proton vertex. The solid line represents data ob- 
tained with a full target and the dashed line represents data 
taken with an empty target. The two inner peaks are events 
produced in the target walls, the two outer peaks represent 
protons produced in the superinsulation of the target cell and 
in its axial heat shield. The range — 7 cm < z < 7 cm (vertical 
lines) has been chosen to select the pp events. 



(where z is measured along the beam line) of the vertex 
more than 7 cm away from the center of the target, as 
shown in Fig. 3, reduced this background to less than 1% 
[19]. 

The particle-detection efficiency of the CLAS is not 
uniform and constant throughout its volume. At the 
edges of the active regions, delimited by the shadows of 
the six superconducting coils, the acceptance decreases 
and varies rapidly. In order to avoid errors, including 
poorly reconstructed tracks in the low-acceptance re- 
gions, a set of fiducial cuts, empirically determined, has 
been applied both to the momenta (pi,P2 > 300 MeV/c, 
p ~ 300 MeV/c being the CLAS detection threshold for 
protons) and to the polar and azimuthal angles (6, <p) of 
the protons. The requirement of having the two protons 
in two different sectors of CLAS has also been applied, in 
order to avoid inefficiencies in the reconstruction of close 
tracks. The angular coverage for the accepted protons is 
shown in the light grey areas of Fig. 4. 

Since the photon energy and the four-momenta of the 
two detected protons are known, and thus the ppn kine- 
matics is completely determined, a missing-mass anal- 
ysis can be performed to identify the neutron. Fig- 
ure 5 shows the missing-mass distribution of the system 
7 3 He — > ppX. The first peak corresponds to the missing 
neutron, the second one to the other competing reaction 
channels, such as those producing pions which had not 
been detected by the CLAS, e.g., 7 3 Hc — > pp(mr°) or 
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FIG. 4: Angular coverage for the identified protons. The gray 
areas represent the fiducial regions of the six CLAS sectors 
inside which the protons for the present analysis have been 
accepted. 



FIG. 6: Examples of missing-mass histograms fitted with a 
Gaussian curve plus an exponential (solid curve) for 0.43 < 
E 1 < 0.45 GeV (al) and 1.13 < £ 7 < 1.15 GeV (a2), for 
0.08 < p n < 0.10 GeV/c and 0.45 < E 1 < 0.55 GeV (bl) and 
0.42 < p n < 0.44 GeV/c and 0.75 < £ 7 < 0.85 GeV (b2), 
and for -0.88 < cos n < -0.84 and 0.35 < £ 7 < 0.45 GeV 
(cl) and 0.72 < cos0„ < 0.76 and 0.95 < £ 7 < 1.05 GeV 
(c2). The background alone is shown as the dashed curves. 



7 3 He — > pp(pTT~). About 25% of the two-proton events, 
~ 5 million events, are thus identified as belonging to the 
ppn channel. 



3500 I 
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FIG. 5: Missing mass of the 7 3 He — » ppX system, for a sub- 
set of the selected pp events. One can easily distinguish the 
peak at the neutron mass at about 0.94 GeV/c 2 (a ~ 0.017 
GeV/c 2 ) from the competing reaction channels. 



The momentum of each detected proton was corrected 
for its loss of energy while passing through the cryogenic 
target material, the target walls, the carbon-fiber scat- 
tering chamber and the start-counter scintillators. 



B. Background Subtraction 

After channel identification, the data were binned in 
photon energy, particle momentum, and particle angle. 
For each of these bins a histogram of the two-proton 
missing-mass distribution was accumulated. Each pp 
missing-mass histogram was fitted with a Gaussian curve 
plus an exponential in order to reproduce the neutron 
peak and the background underneath it. The background 
is due both to misidentified or badly reconstructed pro- 
tons and to the tail from competing reaction channels 
(see Fig. 5). Once the parameters of the fit are extracted, 
the yield is given by the area under the Gaussian curve. 
In this way, the contribution of the background is ex- 
cluded. Some examples of the quality of these fits for 
various bins in photon energy, neutron momentum, and 
neutron angle, chosen to be typical of the character of 
the data for various conditions, are shown in Fig. 6. The 
background-to-signal ratio varies from less than 1% to 
8%, depending on the kinematics. 

To estimate the systematic uncertainty introduced by 
the fitting procedure used to subtract the background 
from the pp missing mass, the yields obtained with two 
kinds of fitting functions for the background (exponential 
and polynomial) have been compared with each other 
[19]. The deviations are, on average, of the order of 2%. 
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C. Efficiency 

Since the neutron is reconstructed using the missing- 
mass technique, the detection efficiency for this chan- 
nel is given by the probability of correctly detecting and 
identifying two protons in the CLAS. This has been eval- 
uated with the aid of a Monte-Carlo simulation. The 
ppn events, generated according to the three-body phase- 
space distribution, were processed by a GEANT-based 
code simulating the response of the CLAS, and were 
reconstructed and analyzed using the same procedure 
adopted for the experimental data. The efficiency in- 
side the CLAS fiducial region for a given kincmatical bin 
At is defined as 



e(Ar) = 



N r 



(1) 



where At lies inside of the CLAS fiducial region, N r is the 
number of reconstructed events within At, and Nq is the 
number of events generated within At. The efficiency so 
computed is more or less constant as a function of photon 
energy, momentum, and angles, and its average value is 
slightly less than 95%. 

In order to evaluate the systematic uncertainties in the 
efficiency for detecting the ppn events in the CLAS, the 
results obtained with the phase-space distribution have 
been compared with the efficiency computed with three 
other event distributions [19]. The result of the calcula- 
tions of the efficiency inside of the CLAS fiducial region 
turns out to be independent of the model used to simulate 
the reaction, apart from the effect of bin migration due 
to the finite resolution of the detector, which has been 
found to be small. The resulting systematic uncertainty 
was determined to be no greater than 5% [19]. 



D. Cross Sections and Normalization 

Three kinds of CLAS-integrated cross sections have 
been measured and are reported here. They are 



total cross sections, defined as 



Le 



(2) 



• semi-differential cross sections with respect to mo- 
mentum, defined as 



da N, 



ppn 



dp ApLe 



(3) 



semi-diffcrcntial cross sections with respect to cos 8, 
defined as 



A, 



da 



dn 27rA(cos0)Le ' 



ppn 



(4) 




FIG. 7: Tagging efficiency as a function of T-counter num- 
ber measured in one particular low-flux run. The average 
efficiency is about 70%. 



where N ppn is the number of events in the bin, e is the 
detection efficiency defined in Section IIIC, and L is the 
luminosity, which is defined as: 



L = N. 



PzNa 
A 



(5) 



where p = 0.0675 g/cm 3 is the density of the target, 
z = 14.0 cm is the effective target length, A is the atomic 
mass of the target (A = 3.016 g/mol), is Avogadro's 
number, and A 7 is the number of incident photons. 

The systematic uncertainties in the target length and 
density arc of the order of 2%. The photon flux was 
measured by integrating the tagger rate over the data- 
acquisition lifetime. The tagging efficiency was measured 
during low-flux runs, using the lead-glass total absorption 
detector. For each T-counter i, the tagging efficiency is 
defined as [21]: 



T eff (z) = (T i .TAC)/Tf 



(6) 



where Tj • T AC is the rate of coincidences between tagger 
and total absorption counter, and T? aw is the rate in the 
tagger alone. A typical tagging efficiency spectrum, as a 
function of T-Counter number, is shown in Fig. 7. To es- 
timate the systematic uncertainty for the photon flux, the 
variations with time of the tagging efficiency and of the 
proton yield normalized to the photon flux for each tag- 
ger scintillation counter have been studied. The resulting 
systematic uncertainty is, on average, approximately 6% 
[19]. The values of the systematic uncertainties in the 
measured cross sections are summarized in Table I. The 
luminosity, integrated over the entire running time and 
over the full photon-energy range, was L ~ 8.7 x 10 35 
cm -2 for this experiment. 



IV. MODEL CALCULATION 

As mentioned in the introduction, the only theoretical 
model currently available for calculation of the cross sec- 
tion for the three-body photodisintegration of 3 Hc in the 
GeV energy region is the one by Laget. In this model, the 
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TABLE I: Systematic uncertainties in the measured cross sec- 
tions. The total is the sum in quadrature of the individual 
uncertainties. 



Quantity 


Uncertainty 


Target length and density 


2% 


Background subtraction 


2% 


Detection efficiency 


5% 


Photon flux 


6% 


Total 


8% 



five-fold differential cross section in the laboratory sys- 
tem for the 7 3 He — > ppn reaction is connected through a 
Jacobian to a reduced cross section 



d 5 <r 



dpdflidfl2 



Expi\E n p\ - E 2 p7i ■ P2I 

d 5 a re d 
(dQi)cmdpndn n 



(7) 



where (Ei,p{), (E 2 ,P2), and {E n ,pX) are, respectively, 
the four-momenta of the two outgoing protons (1 and 2) 
and the neutron in the laboratory frame, and p and Q are 
the proton momentum and the total energy measured in 
the center-of-mass frame of the two protons. 

The reduced cross section depends on the transition 
amplitude T(7 3 i7e — > ppn) [4, 5]: 



d 5 er 



red 



(dQi) cm dp n dQ r 



ppi i 



T|>F 



He I 
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The fully antisymmetrized 3 He bound-state wave func- 
tion \^>3He) is the solution of the Faddeev equations [22] 
for the Paris potential [23]. ft is expanded in a basis 
where two nucleons couple to angular momentum L, spin 
S 1 , and isospin T, the third nucleon moving with angu- 
lar momentum I. Each component is approximated by 
the product of the wave functions, which describe the 
relative motion of the two nucleons inside the pair and 
the motion of the third nucleon [24]. Fermi-motion ef- 
fects are taken fully into account in the two-body ma- 
trix element, and partially [3] in the three-body matrix 
element. However, it has been ascertained that the ef- 
fect of the Fermi motion in the three-body matrix ele- 
ment does not significantly affect the results; therefore, 
it has not been implemented in the version of the model 
which has been used here with the Monte-Carlo proce- 
dure in order to avoid prohibitive computation time. All 
of the S, P, and D components of the 3 He wave function 
arc included. The energy and momentum arc conserved 
at each vertex, and the kinematics is relativistic. The 
continuum final state ppn ) is approximated by a sum 
of three-body plane waves and half-off-shcll amplitudes 
(which are the solutions of the Lippman-Schwinger equa- 
tion for the Paris potential) where two nucleons scatter, 
the third being a spectator. Only S-wave NN scattering 
amplitudes have been retained in the version used in this 
work. The antisymmetry of the continuum final state 




FIG. 8: Diagrams used in Laget's model [3, 8-11] in the cal- 
culation of the 3 He(7,pp)n cross sections: (a) IN absorption 
mechanism, (b) lN+Final State Interactions (FSI), (c) 2N 
absorption, (d), (e) and (f) 2N + FSI, (g) and (h) 3N mecha- 
nisms, and (i) 3N+FSI. The open circles represent full transi- 
tion amplitudes (T matrices); the filled circles are 7 AW and 
ttNN Born terms. 



is achieved by exchanging the roles of the three nucle- 
ons [10]. The transition amplitude T is expanded in a 
truncated series of diagrams that are thought to be dom- 
inant. These diagrams, that were thought to include the 
most likely one-, two- and three-body mechanisms, are 
computed in momentum space. Among all the possible 
three-body mechanisms, meson double scattering is the 
most likely to occur. The Feynman diagrams included 
in the present version of the model are shown in Fig. 8. 
The open circles represent the full transition amplitudes 
(T matrices), which have been calibrated against the cor- 
responding elementary channels, and the filled circles are 
just the 7A/W and irNN Born terms. 

The first two diagrams, (a) and (b), describe one-body 
photoabsorption; (c), (d), (e), and (f) represent two-body 
processes [25, 26]; and (g), (h), and (i) are three-body 
mechanisms, with two-meson (ir or p) exchange. Pion 
absorption by a T = 1 (pn or pp) pair has been found 
experimentally to be strongly suppressed [27], at least 
at low energies, and has not been included in the model 
at this stage. The 3A^ absorption mechanism shown in 
diagram (g) represents the primary 3iV process for the 
3 He(7,pp)n reaction. Above the photon energy corre- 
sponding to the pion-production threshold, the calcula- 
tion does not contain any free parameters, since all of the 
basic matrix elements have been fixed independently us- 
ing relevant reactions induced on the nucleon and on the 
deuteron [4, 5]. The calculated cross section involves a 
logarithmic singularity associated with the on-shcll prop- 
agation of the "first" exchanged pion, which shows up, 
and moves when the photon energy varies, in a well de- 
fined part of the phase space. Below the pion threshold, 
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both exchanged pions are off their mass shells, and the 
three-body exchange currents can be linked by gauge in- 
variance to the corresponding three-body forces [3, 8-11]. 

All model calculations discussed in the following sec- 
tions have been performed with Monte-Carlo sampling 
over the CLAS geometry to produce cross sections that 
can be compared with the experimental results. The 
small-scale structures which are seen in some of the model 
results result from this Monte Carlo treatment, although 
the major structures are real features of the model cal- 
culations. 



RESULTS 



A. Cross Sections Integrated over CLAS 




T P ,/T 



The use of a triangular Dalitz plot is very suitable to 
look for the deviations of an experimental distribution 
from pure phase-space predictions and to identify corre- 
lations between three final-state particles. In particular, 
this technique can be used to identify and select the re- 
gions of the phase space where three-body processes can 
be dominant. If T p i, T p %, and T n are the center-of-mass 
kinetic energies of the two protons and the neutron, re- 
spectively, and T is their sum, we can define the Carte- 
sian coordinates of the triangular Dalitz plot as: 



1 T, 



pi 



V 



V3 

Tn 

T ' 



T 



and 



Figure 9 shows the distribution of the ppn events on the 
Dalitz plot after applying the selection cuts. The wide 
acceptance of the CLAS allows us to fill the physically 
accessible region — delimited by the boundary circle — 
almost completely. 

The shading of the boxes indicate the yield of the ob- 
served ppn events. Areas of increased yield are visible 
where the T p \ and T v i axes intercept the boundary cir- 
cle, as well as where T n w 0. These areas correspond to 
quasi-two-body breakup and neutron-spectator kinemat- 
ics, respectively; they are discussed in detail in Sections 
VA4 and VA2 below. The depletion areas in the upper 
left and upper right sides of the circle correspond to the 
kinematics where one of the protons has low momentum 
(Pi>P2 < 300 MeV/c) and therefore is not detected by the 
CLAS. The central top area where the two protons are 
emitted in nearly the same direction is excluded by the 
requirement of detecting the two protons in two differ- 
ent sectors (see Section III A above) . The central region, 
near the intersection of the three axes, consists of events 
where all three nucleons have nearly equal energies, and 
is called the "star" region (see Section V A 3 below) . 

In the following sections, CLAS-integrated cross sec- 
tions for the full acceptance and for the three selected 
kinematics listed above, each chosen to illustrate its two- 



FIG. 9: Triangular Dalitz plot for the ppn data. T p i, T P 2, 
and T n are the center-of-mass kinetic energies of the three 
nucleons. 



body or three-body character, are presented and com- 
pared with distributions obtained both with three-body 
phase space and with the results of the Laget model. 



1. Full CLAS acceptance 

The ppn total cross section integrated over the CLAS 
acceptance has been measured as a function of the inci- 
dent photon energy E y . The photon-energy spectrum, 
ranging from 0.35 GeV to 1.55 GeV, has been divided 
into 60 bins, each 0.02 GeV wide. The results are shown 
in Fig. 10. The cross section, ranging between 10 /xb 
and 0.01 /xb, decreases almost exponentially as the pho- 
ton energy increases. Fitting the data with an expo- 
nential function cr(£ , 7 ) cx e~ bE ~< yields a slope b ~ 5.3 
GeV -1 . The data are compared with the full calculation 
(solid curve), including one-, two-, and three-body mech- 
anisms, as well as with the results for the one- and two- 
body mechanisms only (dashed curve), and the three- 
body mechanisms only (dotted curve), as shown in Fig. 8. 
It is important to note that the theoretical curves repre- 
sent absolute cross sections calculated within the CLAS 
acceptance - they are not normalized to the data. The 
results of the model calculations that do not include the 
three-body mechanisms are almost a factor of ten smaller 
than the data at lower energies, while they approach the 
data as the photon energy increases and exceed the data 
at higher energies. The full-model results agree better 
with the data, but still are too low at low energies and 
too high at high energies. 

Figure 11 shows the partial differential cross section 
as a function of neutron momentum p n , for twelve 0.1- 
GeV-wide photon energy bins. The data are compared 
with phase-space-generated event distributions (dotted 
curves) normalized in each energy bin in order to match 
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FIG. 10: Total ppn cross section integrated over the CLAS 
acceptance plotted as a function of photon energy on a log- 
arithmic scale for the full E y range. The ppn cross section 
(circles) is compared with Laget's full model (solid curve), 
with the model result without the three-body mechanisms 
(dashed curve), and with the one including only three-body 
processes (dotted curve). The error bars include statistical 
and systematic uncertainties, as in all the following experi- 
mental distributions. 



the area under the experimental distribution, with the 
results of Laget's full model (solid curve), and with the 
model with no three-body mechanisms included (dashed 
curve). The neutron momentum distributions are related 
to the projection of the data in Fig. 9 onto the T n axis. 

In the photon-energy range between 0.35 and 0.95 
GeV, the data show a broad central distribution in the 
middle of the neutron momentum spectrum {e.g., at 
about 400 McV/c for £ 7 = 0.4 GeV and 500 MeV/c 
for E 1 > 0.5 GeV), which is reproduced reasonably well 
by the phase-space distribution (better at low photon 
energies than at high energies). Up to about 0.6 GeV 
a comparison of the data with the shape of the model 
results reveals the presence of three-body mechanisms. 
In the middle range of neutron momentum, two-body 
mechanisms are seen to contribute increasingly starting 
from E 1 = 0.65 GeV. These contributions stem from 
low-energy S-wave np rescattering, which causes the in- 
creased yield in the quasi-two-body kinematics, corre- 
sponding to the the left and right sides of the Dalitz plot 
(Fig. 9). This yield projects onto the middle range of the 
neutron-momentum distribution. 

A peak, roughly 0.04-GcV/c wide, is observed at a 



neutron momentum of about 0.12 GeV/c, independent 
of the photon energy. The relative strength of this peak 
increases with increasing photon energy, but it is not ac- 
counted for by the three-body phase-space distribution. 
However, this structure is expected by the model, and it 
is predicted to be largely due to two-body mechanisms. It 
reflects the Fermi distribution of the spectator neutron. 
This feature has been exploited to select the neutron- 
spectator kinematic region, as is explained in Section 
VA2 below. 

At photon energies from about 0.9 to 1.2 GeV and high 
neutron momenta, a third structure appears in the data, 
which is present neither in the phase-space distribution 
nor in the (l+2)-body model results, but is predicted by 
the full model. This structure can therefore be considered 
to be a signature of three-body mechanisms as well. 

The differential cross section as a function of the co- 
sine of the neutron polar angle cosf5„ in the lab system 
is plotted in Fig. 12, for twelve 0.1-GeV-wide photon- 
energy bins, between 0.35 GeV and 1.55 GeV. The dis- 
tributions are forward-peaked at low-to-intermediate en- 
ergies, while they become flatter for higher E 1 . Their 
shapes are reasonably well reproduced by both phase- 
space and the full- model calculations. 



2. Spectator neutron 

Guided by Fig. 11, the events where the neutron is a 
spectator in the photobrcakup of a proton pair have been 
selected by requiring the condition p n < 250 McV/c. 
These are all the events in the lower neutron-momentum 
peak (within 3cr from its center). 

Figure 13 shows the cross section as a function of pho- 
ton energy integrated over the CLAS for the events sat- 
isfying this condition, compared with the predictions of 
the model. After an initial steep drop, the cross sec- 
tion has an exponential dependence on the photon en- 
ergy above 0.6 GeV, this time with a slope b ~ 4 GeV -1 . 
The agreement between the experimental cross section 
and the model prediction is good only for low energies, 
below 600 MeV. The cross section is clearly driven by 
two-body mechanisms, as expected. 

The differential cross section as a function of 
(cos 9 n ) i a b, which is plotted in Fig. 14 for eight photon- 
energy bins, shows a generally flat distribution. This is 
expected, because in the neutron-spectator kinematics 
the two-body part of the reduced differential cross sec- 
tion is proportional to the neutron-momentum distribu- 
tion p(n) times the center-of-mass differential cross sec- 
tion for the pp-pair breakup [8] : 

= (1 + (3 n cos 6 n )p(p n )-^—(jpp ^ pp). (9) 

/ ^^cmO y Pn Ctil cr n 

Both the (l+2)-body part and the full-model results 
agree fairly well, in shape and magnitude, with the ex- 
perimental distributions up to 600 MeV. At higher en- 
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FIG. 11: Differential cross sections integrated over the CLAS as a function of the neutron momentum in the laboratory frame for 
twelve 0.1-GeV-wide photon-energy bins between 0.35 GeV and 1.55 GeV. The points represent our CLAS data. The error bars 
include both statistical and systematic uncertainties. The dotted curves are the distributions for phase-space events generated 
within the CLAS acceptance and normalized in each energy bin to match the total area of each experimental distribution. The 
solid curves represent the full Laget-model results, while the dashed lines represent the model including one- and two-body 
mechanisms only. For E 1 > 0.95 GeV, the model predictions at p„ < 250 MeV/c (to the left of the vertical dotted-dashed line) 
are scaled by a factor 0.1 to fit in the plots. 
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FIG. 12: Differential cross sections integrated over the CLAS as a function of the cosine of the neutron polar angle in the 
laboratory frame for twelve 0.1-GeV-wide photon-energy bins between 0.35 GeV and 1.55 GeV. The points represent our CLAS 
data. The error bars include statistical and systematic uncertainties. The dotted curves are the distributions for phase-space 
events generated within the CLAS acceptance and normalized in each energy bin to match the total area of each experimental 
distribution. The solid curves represent the full model results, while the dashed lines represent the model including one- and 
two-body mechanisms only. 
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FIG. 13: Cross section integrated over the CLAS for the 
neutron-spectator kinematics plotted as a function of photon 
energy. The data are compared with the predictions of the 
full model (solid curve), the (l+2)-body part (dashed curve), 
the three-body part (dotted curve), and the 1-body part alone 
(dashed-dotted curve). 



crgies, the calculation predicts the contribution of two- 
body mechanisms to be much too large. 

In the neutron-spectator kinematics, the primary 
physics is contained in the angular distribution of the 
1PP ~ * PP subchannel. Figure 15 compares this angular 
distribution with the model. While the magnitude of the 
experimental cross section is well reproduced at low en- 
ergy by the model, the shape of the angular distribution 
is markedly different. The model curve exhibits a min- 
imum at 90°, where the measured differential cross sec- 
tion has a broad maximum. It can be seen from Fig. 13 
that three-body diagrams do not contribute significantly 
to the total cross section, but their interference with the 
two-body diagrams brings the shape of the angular dis- 
tributions closer to the experimental ones. However, this 
effect is not strong enough to cancel the huge contribu- 
tion of the two-body part at high energy. 

Since the pp pair that absorbs the photon has no 
dipolc moment for the photon to couple with, charged- 
meson exchange currents and intermediate-A production 
(Fig. 8, diagrams (c) and (d)) are strongly suppressed 
and one-body mechanisms (diagrams (a) and (b)) and 
related FSI (diagrams (e) and (f)) contribute more signif- 
icantly to the two-body photodisintcgration cross section 

ilPP ~ * PP) ■ At low energy, the one-body amplitude is 
driven by dipole photon absorption, which is suppressed. 



At high energy, it involves all other multipoles and, as 
a result, the corresponding cross section remains almost 
constant. This process probes the relative pp wave func- 
tion at a momentum which increases with the incoming 
photon energy — typically 400 MeV/c at E y = 400 MeV, 
increasing to 1.5 GeV/c at E 7 = 1.2 GeV. Above ~ 0.8 
GeV the pp wave function is not under control, and we 
are reaching the limits of the model, as in the jd — > pn re- 
action [28] . We may have entered a region where quarks 
become the relevant degrees of freedom [28, 29], or per- 
haps a description in terms of Rcggc-typc calculations 
[30] is more suitable. 



3. Star configuration 

The center of the Dalitz triangle corresponds to the 
three particles having equal kinetic energies and their 
three-momentum vectors forming angles of 120° with 
each other (in the ppn center-of-mass frame). For this 
reason, this kinematical arrangement, shown schemati- 
cally in Fig. 16, has been called the star configuration. In 
this region, the three-body mechanisms are expected to 
be dominant because if the momentum is equally shared 
between the three nucleons, the contribution from two- 
body mechanisms is minimized. This is therefore consid- 
ered to be a good place to study three-body mechanisms. 

The events for this kinematics have been selected by 
requiring that the three nucleons satisfy the condition 

|% - 120°| < Sg (10) 

where 

%=arccosf^i^ (11) 
V PiPj ) 

is the angle between the momenta of nucleons i and j, 
in the center-of-mass frame, and the angle Sg, which ex- 
presses the allowed deviation from the pure "star" kine- 
matics, has been chosen to be 15°, as shown in Fig. 17. 

In Fig. 18, the cross section integrated over the CLAS 
for the star configuration is plotted as a function of 
photon energy. It decreases exponentially, with slope 
b ~ 5.8 GeV -1 , as the photon energy increases, much 
more steeply than for the neutron-spectator kinematics. 

As expected from the kinematics, for the star configu- 
ration the contribution of two-body mechanisms is neg- 
ligible, while the bulk of the cross section comes from 
three-body mechanisms. At low energy the model misses 
the experimental cross section by approximately a factor 
of four. The probable reason for this discrepancy is that 
only the Born term and the A-formation term [31] have 
been retained in the calculation of the 7^ — > ttN vertex 
(the upper blob in Fig. 8 (g) and (h)). The addition of 
the contributions of the iV(1520)D 13 , AT(1440)Pu, and 
7V(1535)Sn resonances also might improve the agree- 
ment with the data. At high E 1 , the Blomqvist-Laget 



13 



tr 
40 

— 30 - 

I 20 - 
■a 

10 - 



E = (0.35-0.45) GeV 




: (0.45-0.55) GeV 




: (0.55-0.65) GeV 





a 

To 



a 

"O 
"O 



a 



e 

"O 



- a 



-0.5 





cose n LAB 



0.5 



20 



15 



10 




12 

10 




10 

8 

6 

4 

2 


4 



E = (0.75-0.85) GeV 



E = (0.85-0.95) GeV 



E T = (0.95-1.05) GeV 



E y = (1.05-1.3) GeV 



-0.5 

cos6„ LAB 



0.5 



FIG. 14: Differential cross sections integrated over the CLAS for the neutron-spectator kinematics with respect to cos 9 of the 
neutron in the laboratory frame for photon energies between 0.35 and 1.30 GeV. The data are compared with the results of 
the full model (solid curves) and those of the one- plus two-body-only model (dashed curves), for 0.35 < < 0.75 GeV only, 
because at higher energies the model calculations differ by more than an order of magnitude from the data. 
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FIG. 15: Differential cross sections integrated over the CLAS for the neutron-spectator kinematics with respect to cos6 pp of 
the proton in the pp-pair center-of-mass frame for photon energies between 0.35 and 1.30 GeV. The data are compared with the 
results of the full model (solid curves), and those of the one- plus two-body-only model (dashed curves), for 0.35 < E-y < 0.75 
GeV, because at higher energies the model calculations differ by more than an order of magnitude from the data. 
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FIG. 16: Kinematics of the star configuration in the ppn 
center-of-mass frame. The angles 9* , between the normal vec- 
tor to the star plane and the photon-beam direction, and <f>* , 
the neutron azimuthal angle in the star plane, define the re- 
action. 
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FIG. 17: Dalitz plot for the CLAS ppn events selected for the 
star configuration. 

Born term matches the Regge amplitudes [32] that re- 
produce the 7-/V — > ttN cross section in this energy re- 
gion. The pion-rescattering amplitude (Fig. 8 (g)-(i)) is 
parametrized in terms of partial waves up to and includ- 
ing G- waves. 

The differential cross section as a function of cos 9* , the 
cosine of the angle between the incident photon and the 
normal vector to the three-nucleon center-of-mass plane 
(see Fig. 16), is plotted in Fig. 19 for eight photon-energy 
bins between 0.35 GeV and 1.30 GeV. Since the two out- 
going protons are indistinguishable, the orientation of the 
normal vector to the star plane, pi x p 2 , is arbitrary. 
Thus, the distribution is symmetric around cos#* = 0. 
The shape of the cross section is very well reproduced by 
phase space at low energy, while at high energy the model 
better reproduces the curvature of the experimental dis- 
tribution. At all energies the three-body mechanisms are 
dominant. 

Figure 20 shows, for eight photon-energy bins between 
0.35 GeV and 1.30 GeV, the differential cross section as 
a function of the angle <fi* between the neutron direction 
in the star plane and the projection of the photon-beam 
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FIG. 18: Cross section integrated over the CLAS acceptance 
for the star configuration plotted as a function of incident 
photon energy. The CLAS data are compared with the pre- 
dictions of the full model (solid curve), to the one- plus two- 
body-only part (dashed curve), and to the three-body-only 
part of the model (dotted curve). 

direction in the same plane (see Fig. 16). As is the case 
for 8* , the angular distribution is symmetric, here around 
180°. It also follows a phase-space distribution, except 
for E 1 > 0.95 GeV, and its shape (but not its magnitude) 
is reproduced fairly well by the model as well. Again, 
three-body mechanisms are seen to be dominant. 

The photoproduced pion described by the diagrams (g) 
and (h) of Fig. 8 can propagate on-shell, since the avail- 
able energy is larger than the sum of the masses of the 
pion and the three nuclcons. This causes the development 
of a logarithmic singularity in the thrcc-nuclcon ampli- 
tude, which should enhance the contribution of three- 
body mechanisms. The effect of this singularity can be 
seen in Fig. 21, in which is plotted the cross section dif- 
ferential in m 2 x /m 2 , where mx, defined from the relation 

m x = {E 1 +m p -E n f -(k^-p^f, (12) 

is the missing mass in the 7p — > n + n reaction, assuming 
that the proton is at rest. 

At photon energies above about 0.6 GeV, the pion sin- 
gularity appears clearly (m x /m 2 ~ 1) in both the ex- 
perimental distributions and the model results. At high 
energy, the magnitudes of the two peaks are compara- 
ble but the shift of the theoretical one with respect to 
the experimental one reflects the approximate treatment 
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FIG. 19: CLAS-integrated differential cross sections with respect to cos#* for the star configuration. The data, for photon 
energies between 0.35 GeV and 1.30 GeV, are compared with the full-model results (solid curves) and the one- plus two-body- 
only part (dashed curves). The dotted curves are the phase-space distributions multiplied, for each photon-energy bin, by the 
constants used to normalize the full-Dalitz cross sections. 
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energy bins. The dotted lines represent the phase-space pre- 
dictions, multiplied by the constants used to normalize the 
full-Dalitz cross sections, while the solid curves are the full- 
model results. 



of Fermi motion effects in the model. At lower energy, 
the theoretical peak is smaller than in the experiment. 
The inclusion of higher-lying resonances in the sequen- 
tial scattering amplitude in the model will enhance the 
peak near m\jm\ — 1, but will probably not fill the gap 



around 



'A" 



/mi 



-15 for £L = 400 MeV. 



These findings indicate a deviation from the sequential 
rescattcring three-body mechanisms, which may be a hint 
in the search for more genuine three-body processes. 



4- Quasi-two-body breakup 

The third region of the Dalitz plot examined corre- 
sponds to the quasi-two-body breakup, where a proton 
and an unbound deuteron (a pn pair) are emitted back- 




FIG. 22: Dalitz plot for our CLAS ppn events selected as 
quasi-two-body breakup. 



to-back in the center-of-mass frame. For this kind of 
event, one of the two protons (pi) is emitted with 2/3 
of the total available energy, and the pn pair travels in 
the opposite direction, with 1/3 of the total energy, and 
with T P 2 = T n = gT. This kinematics corresponds to 
the events in the two populated areas shown in Fig. 22. 
These areas have been selected by requiring that the an- 
gle between the high-energy proton and each of the other 
two nucleons be close to 180°, and that the difference 
between the energies of the two low-energy nucleons be 
small. Using the formalism defined above, 



Wj-180°| < 20°, 
|0 P i„-18O°| < 20°, 



(13) 



and 



IT, 



p2 



T 



< 0.15 



for the events on the right side of the Dalitz plot (where 
the proton labeled pi has higher energy), and 



\6p2n 



180° | < 20° 
180° I < 20° 



(14) 



and 



T — T 
- L pl ± n 



< 0.15 



for the events on the left side of the Dalitz plot. Since 
protons "1" and "2" arc indistinguishable, the two re- 
gions of the Dalitz plot are equivalent. 

In Fig. 23, the CLAS-integrated cross section for this 
process is plotted as a function of photon energy. It de- 
creases exponentially with a much steeper slope than for 
the other kinematics (b ~ 7.3 GeV" 1 ). The full- model 
result is in good agreement with the experimental cross 
section only for the low part of the photon energy range, 
and seriously underestimates it above about E 1 = 0.55 
GeV. The (l+2)-body calculation gives a cross section 
that is smaller than the data by a factor of five or more 
for all photon energies. However, this kinematic region 
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FIG. 23: Cross section integrated over the CLAS for the 
quasi-two-body breakup plotted as a function of photon en- 
ergy. The data are compared with the predictions of the 
full model (solid curve), the (l+2)-body calculation (dashed 
curve), and the three-body-only calculation (dotted curve). 
The full-model calculation agrees quantitatively with our ex- 
perimental results only up to about 0.55 GeV. 



is expected to be strongly influenced by final-state in- 
teractions (FSI) [14]. Only S-wave NN scattering has 
been included in the model calculation. Furthermore, a 
factorization approximation has been made to estimate 
the nine- fold integral in Fig. 8, graph (i). A full treat- 
ment, in the terms of Ref. [33], might help to reduce the 
discrepancy between the data and the model predictions. 

It also turns out that the logarithmic singularity in 
the two-step sequential scattering (Fig. 8 (g) and (h)) 
moves in the Dalitz plot as the photon energy varies. At 
lower photon energies, around E~ = 500 McV, it coin- 
cides with the part of the Dalitz plot where the quasi- 
two-body events are located and where the amplitude 
includes a significant contribution from FSI as well. As 
the photon energy increases, the singularity moves to- 
wards the top of the Dalitz plot, and the contribution 
of sequential scattering to the quasi-two-body cross sec- 
tion becomes negligible. Here, the difference between the 
experimental cross section and the full-model result is a 
strong hint of a possible contribution of other three-body 
mechanisms that do not reduce to sequential scattering. 

In Fig. 24, the differential cross section is plotted as 
a function of the cosine of the polar angle of the higher- 
energy proton in the three-body center-of-mass frame. 
Data from eight photon-energy bins between 0.35 and 



1.30 GeV are shown. The experimental cross section 
shows a forward peak whose relative strength grows with 
increasing photon energy. This feature is also seen in 
the (l+2)-body model and in the full calculation for 
Ey > 0.55 GeV. The predicted strength of the forward 
peak is, however, much too small to match the data. For 
lower energies, the full calculation predicts a cross-section 
enhancement at backward angles that is not seen in the 
data. 



B. The ppn "Three-Body" Cross Section 

Previous experiments measuring the 7 3 He — ► ppn 
channel in an extended part of the phase space have been 
performed with the DAPHNE [18] and TAGX [16] detec- 
tors. Except for differences in the cj> coverage, the CLAS 
event-selection cuts are very similar to the other two ex- 
periments, as seen in Table II; however, differences in 
the selection criteria of the three-body events exist be- 
tween the TAGX experiment on the one hand and the 
DAPHNE and CLAS experiments on the other. 

The ppn "three-body" cross section is defined as: 



Cr3body{E^) 



Nobody (Ej) 



N y (E y ) Acc(E^) 



(15) 



where Nobody is the number of events extracted by ap- 
plying the selection cuts given in Table II and Acc is the 
acceptance of the CLAS detector for the ppn events calcu- 
lated with the phase-space Monte-Carlo simulation. The 
low-momentum neutrons (p„ < 150 MeV/c) have been 
excluded in order to select only those events for which 
all three particles participate in the reaction, thus di- 
minishing the importance of two-body processes [16, 18]. 
In this kinematics, the phase-space result describes the 
process reasonably well. 

Figure 25 shows a^body as a function of the photon en- 
ergy E 1 . The full circles represent our CLAS data, the 
empty triangles the data of the TAGX Collaboration [16], 
and the empty squares the results obtained in the exper- 
iment carried out at MAMI with the DAPHNE detector 
[18]. The error bars on the CLAS data are statistical 
only. The systematic uncertainties delineated in the pre- 
vious section are shown by the vertical lines in the upper 
part of the figure. 

In the overlap region of the three experiments from 
0.35 to 0.80 GeV, the CLAS data are in good agreement 
with the DAPHNE results, but differ from the TAGX 
cross sections by about 15%, most likely due to the above- 
mentioned difference in the three-body event selection. 
Above 0.80 GeV, no previous data are available. 

The phase-space extrapolation to the unmeasured re- 
gions has been done only for comparison with the previ- 
ous experiments, which adopted the same procedure to 
extract cr 3 body 
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FIG. 24: Differential cross sections integrated over the CLAS for the quasi-two-body breakup with respect to cos# of the 
high-energy proton in the center-of-mass frame for photon energies between 0.35 and 1.30 GeV. Our data, for 0.35 < E 1 < 0.75 
GeV, are compared with the results of the full model (solid curves) and of the (l+2)-body-only model (dashed curves). 
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TABLE II: Selection cuts applied to the TAGX, DAPHNE, and CLAS 7 3 He — » ppn experiments in order to extract the 
"three-body" total cross section. 



TAGX 


DAPHNE 


CLAS 


15° < 9 P1 . P2 < 165° 
0° < <t> PUPa < 40° 
Ppi,p2 > 300 MeV/c 
"Non-spectator" neutron 


22° < 9 P1 , P2 < 158° 
0° < <P P1 , P2 < 360° 
Ppi.pi > 300 MeV/c 
p n > 150 MeV/c 


15° < 6 pltP2 < 125° 
CLAS <f> fiducial cuts 
P P1 ,P2 > 300 MeV/c 
p„ > 150 MeV/c 
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FIG. 25: Total "three-body" cross section as defined by 
Eq. (15) for the 7 3 He — > ppn reaction plotted as a function of 
photon energy. The CLAS data (full circles) are compared 
with the results from DAPHNE [18] (empty squares) and 
TAGX [16] (empty triangles). The error bars on the CLAS 
experimental points are statistical only. The CLAS system- 
atic uncertainties are represented by the vertical bars in the 
upper part of the figure. 



VI. SUMMARY AND CONCLUSIONS 

The three-body photodisintcgration of 3 Hc has been 
measured with the tagged-photon beam and the CE- 
BAF Large Acceptance Spectrometer in Hall B at the 
Thomas Jefferson National Accelerator Facility in the 
photon-energy range between 0.35 and 1.55 GeV. This 
measurement constitutes a wide-ranging survey of two- 
and three-body processes in the 7 3 He — ► ppn reaction 
channel, as a consequence of the high statistics and large 
kinematic coverage obtained with the CLAS. 

Total and partially integrated differential cross sections 
for the full ppn data set and for selected kinematics were 
extracted and are compared with phase-space distribu- 



tions and with the predictions of the diagrammatic model 
of Laget. This model reproduces some of the main trends 
of the experimental energy distributions, and for these 
cases can be taken as a qualitative guide to understand- 
ing the reaction mechanisms. 

From the analysis of the neutron-momentum distribu- 
tion for the full Dalitz plot, the kinematic region corre- 
sponding to the photodisintegration of a pp pair in the 
presence of a spectator neutron has been identified. Here, 
the effects of two-body absorption mechanisms dominate 
and the model results are very close to experiment at 
low energy, up to E 1 = 600 MeV. At higher energies, 
the discrepancy, which increases with energy, might be a 
hint that we are approaching the limit of models based 
on meson and baryon degrees of freedom. 

A strong contribution of three-body sequential meson- 
absorption mechanisms is manifested over all the avail- 
able phase space, but most especially in the star kine- 
matics, the spatially symmetric configuration of the three 
final-state nucleons. These events are dominated by the 
coupling to the A resonance, and they strongly confirm 
its role in three-body forces. The deviations from the 
predictions of the diagrammatic model point not only 
toward the necessity of implementing processes which in- 
volve higher-lying baryonic resonances, but also toward 
possible additional three-body mechanisms beyond se- 
quential scattering. 

The 47r-integrated "three-body" cross section is in 
excellent agreement with previous experimental results 
from DAPHNE up to 800 MeV. For the first time we 
now have provided access to a higher energy range, up to 
1.5 GeV. 

This work breaks new ground in the experimental 
study of the three-body photodisintegration of 3 Hc. 
However, before making contact with the elusive three- 
body forces, it calls for a more complete treatment of 
three-body mechanisms which go beyond the dominant 
sequential meson exchange and A formation in the in- 
termediate energy range, and which take into account 
possible coupling with partonic degrees of freedom in the 
highest energy range. 
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